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ABSTRACT 

We explore the cosmic evolution of radio AGN with low radio powers (Li 4qhz J$ 5 x 10 25 W Hz -1 ) 
out to z — 1.3 using to-date the largest sample of ~ 600 low luminosity radio AGN at intermediate 
redshift drawn from the VLA-COSMOS survey. We derive the radio luminosity function for these 
AGN, and its evolution with cosmic time assuming two extreme cases: i) pure luminosity and ii) pure 
density evolution. The former and latter yield L* oc (l + z)°- 8±01 , and <&» oc (l + z) 11±01 , respectively, 
both implying a fairly modest change in properties of low radio-power AGN since z — 1.3. We show 
that this is in stark contrast with the evolution of powerful (Li.4GHz > 5 x 10 25 W Hz -1 ) radio AGN 
over the same cosm ic time interval, constrained using the 3CRR, 6CE, and 7CRS radio surveys by 
IWillott et ail (|2001[ ). We demonstrate that this can be explained through differences in black hole 
fueling and triggering mechanisms, and a dichotomy in host galaxy properties of weak and powerful 
AGN. Our findings suggest that high and low radio-power AGN activity is triggered in different stages 
during the formation of massive red galaxies. We show that weak radio AGN occur in the most massive 
galaxies already at z ~ 1, and they may significantly contribute to the heating of their surrounding 
medium and thus inhibit gas accretion onto their host galaxies, as recently suggested for the 'radio 
mode' in cosmological models. 

Subject headings: galaxies: fundamental parameters - galaxies: active, evolution - cosmology: obser- 
vations - radio continuum: galaxies 



1. INTRODUCTION 

1.1. AGN feedback: The impact of radio sources on 
galaxy formation 

Radio activity from active galactic nuclei (AGN) has 
recently been invoked in cosmological models as a signifi- 
cant ingredient in the process of galaxy formation ('AGN 
feedback'). Given that, in the past, cosmological models 
led to a systematic over-prediction of the high-mass end 
of the galaxy stellar mass function (e.g. IWhite fc Reesl 
119781 IWhite fc Frenklll991l ). the implementation of gas 
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heating through energetic radio outflows - the 'radio' 
mode 16 - managed to reproduce well many observed 
galaxy properties (e.g. the galaxies' stellar mass func- 
tion, color bimodality; Croton et al. 2006, Bower et al. 
2006). The particular choice of a 'radio mode' as the rel- 
evant heating mechanism has been motivated by many 
observational results verifying the interplay between the 
emission of radio galaxies in galaxy clusters, and the clus- 
ter's hot X-ray emitting gas on large-scales (e.g. Fabian 
et al. 2003, Forman et al. 2005). 

In the centers of galaxy clusters the radiative cooling 
time of the intra-cluster medium (ICM) is shorter than 
the age of the cluster. Thus, the central ICM is expected 
to be significantly colder compared to outskirt-regions 
(such clusters are referred to as 'cooling flow clusters'). 
However, generally the expected cool X-ray phases are 
not observed, and this is known as the 'cooling flow prob- 
lem'. The most likely solution to this problem is thought 
to be radio galaxies as it is usually observed that the syn- 
chrotron plasma ejected by their radio jets inflates bub- 
bles in the h ot X-ray emi tting cluster gas and thereby 
heats it (see lFabianlTl994l for a review) . In theoretical 
models, a similar interplay on smaller scales has been 
assumed to be at work between the radio outflows of a 
galaxy and its surrounding hot gas halo. 

The first observational evidence for the 'radio mode' 
feedback in the cont e xt of galaxy formation has been 
found by iBest et ail (j2006) . Based on X-ray and ra- 
dio obse rvations of a loc al sample of massive elliptical 
galaxies (jBest et al.ll2005j ) they have shown that the ra- 
dio source heating may indeed balance the radiative en- 
ergy losses from the hot gas surrounding the galaxy, as 
postulated in the models. Based on an independent 
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sample of radio galaxies in galaxy clusters (studied by 
iBzrzan et aLll2004f) iBest et al.l (|2006l ) have correlated the 
mechanical heating provided by the radio sources on their 
surrounding medium with monochromatic 20 cm radio 
power. Combini ng this with the local radio AGN lumi- 
nosity function ([Best et al.l I2Q05T ) . they have found that 
the volume-averaged mechanical energy heating rate of 
local radio luminous AGN is about a factor of 10-20 lower 
than predicted by the Croton et al. (2006) model. They 
assigned this difference to a systematic over-prediction 
of the heating rate in the model (by a factor of 2-3) and 
the intention of the model to jointly describe both the 
large- and small- scale heating from radio sources (i.e. 
on both cluster and galaxy scales). More recent calcula- 
tions of the volum e averaged mechanical energy heating 
rate due to AGN (|K5rding et al.ll2008t IMerloni fc Heind 
2008) yield higher v alues at z — compared to the results 
of Best et all (|2006f). and thus closer to the lCroton et al.l 
(20065 prediction. 

To date no clear picture exists on how the AGN ra- 
dio mode feedback works. More observations are needed 
to provide better constraints on the physics of this pro- 
cess, as well as its evolution with cosmic time. Here 
wc attempt to shed light on the latter utilizing a large 
statistical sample of radio AGN, drawn from the VLA- 
COSMOS survey (jSchinnerer et al.ll200l . 

1.2. Radio luminous AGN and their evolution 

iFanaroff fc Rilevl (|1974 FR hereafter) were the first 
to note that the radio luminous AGN population con- 
sists of two apparently distinct types of sources, re- 
ferred to as FR type I and II, with prominent differ- 
ences in both radio morphology and luminosity. The 
radio emission from FR I radio galaxies is core dom- 
inated, while FR lis are edge-brightened with highly 
collimated large-scale jets. FR lis are typically more 
powerful in the radio than FR Is, and a dividing lumi- 
nosity of Li78MHz ~ 2.5 x 10 26 W Hz" 1 (correspond- 
ing to £i.4GHz ~ 6 x 10 25 W Hz -1 ) 17 has been sug- 
gested (IFa naroff & Rilcvl ll97l ). It has been later shown 
(jLedlow fc Owenlll996h that the FR I/FR II division is 
also a function of the host galaxy optical luminosity (with 
a higher dividing radio luminosity for higher optical host 
luminosity) . 

An alternative way of classifying radio AGN is based 
on the existence of high excitation (HE) emission lines in 
the optical spectra of th eir host galaxies (jHine fc Longairl 
Il979t lLaing et aTH l994). In this scheme, objects with- 
out high-excitation emission lines are referred to as low- 
excitation (LE) radio galaxies, and they are most com- 
mon at low radio luminosities. Almost all FR I radio 
galaxies are LE sources, while optical hosts of the most 
powerful radio sources, i.e. FR lis, usually have strong 
emission lines. It is noteworthy, however, that the cor- 
respondence between the FR class and the presence of 
emission lines is not one-to-one. Many FR II galaxies 
have been found t o be low-excitation radio galaxies (e.g. 
lEvans et al.ll2006h . 

Recently, strong e viden ce (|Evans et al.l 120061 : 
lHardcastle et al.l l2006l l2007ft has emerged support- 
ing the idea that low- and high- excitation radio AGN 

17 A spectral index of a = 0.7 (F v oc u~ a ; where F v is the radio 
flux density) is assumed throughout the paper. 



(hereafter LERAGN and HERAGN, respectively) 
reflect different modes of black hole (BH) accretion. 
Independent studies have shown that i) LERAGN are 
a class of radi o luminous AGN that accrete radiativel y 
inefficiently (jEvans et all 120061 : lHardcastle et al"1l2006h . 
and ii) Bondi accretion of the hot X-ray emitting inter- 
galactic medium (IGM) is sufficient to power the jets of 
low-power radio galaxies in the centers of galaxy clusters 
(based on a sample of nearby galaxies; Allen et al. 
2006). Based on these findings Hardcastle et al. (2007) 
have developed a model in which all low-excitation radio 
sources are powered by radiatively inefficient accretion 
(i.e. at sub-Eddington accretion rates) of the hot phase 
of the IGM, while high-excitation sources are powered 
by radiatively efficient accretion (at ~ Eddington rates) 
of cold gas (that is in general u nrelated to the hot IGM; 
see also IMerloni fc Heinzll2008l) . This model successfully 
explains a variety of properties of radio luminous AGN, 
such as their environments, signs of galaxy mergers in 
the hosts of powerful (high-excitation) radio sources, 
and the break of the radio AGN lum inosity function 
(for details see lHardcastle et al.l 120071 and references 
therein) . 

In the past two decades it has become clear that radio 
luminous AGN evolve differentially: Low-power sources 
evolve less strongly than high-power sources. Numer- 
ous studies of high luminosity radio AGN (L2.7GHZ ^ 

10 25 W Hz -1 sr" 1 corresponding to L i 4qhz ~ 2 x 

10 26 W Hz" 1 ; iDunlop fc Peacock! Il990l : IWillott et all 
l2001h have found a strong positive density evolution 
with redshift of these sources out to a redshift peak of 
z ~ 2. Beyond this redshift their comoving volume den- 
sity starts declinin g (a possible sharp densit y cut-off has 
been suggested bv lDunlop fc Peacocklll990L but has not 
been confirmed by IWillott et al.l l200ir T Such a decline 
would be consistent with the results obtained via optical 
surve ys (ISchmidt et al.lll995h and, recently. X-rav sur- 
veys (jSilverman et al.ll2008l: iBrusa et al.ll2008l ). 

Analyzing the evolution of lower-power radio AGN 
(Li.4GHz > 2 x 10 25 W Hz- 1 ) iWaddington et all (l200lh 
have found a significantly slower evolution of this popu- 
lation, with the comoving volume density turn-over oc- 
curring at a lower redshift (z ~ 1 — 1.5), compared 
to the high-power radio population. Different radio- 
optical surveys are still som ewhat controversial . While 
iClewlev fc Jaxvisl (l200l and lSadler et~aT1 (pOOl find no 
evidence for partic ularly strong evolution out to z ~ 0.7, 
ICowie et al.l (|2004[ ) have found a strong density evolution 
out to z ~ 1.5. Thus, although on average a weaker evo- 
lution of low-power (compared to high-power) radio AGN 
is implied, it is still not very well understood how the low- 
luminosity radio AGN evolve with redshift. In this work 
we use the VLA-COSMOS AGN sample of low luminos- 
ity (96% have Li. 4G Hz < 10 25 W Hz -1 ) radio sources in 
order to comprehensively constrain the evolution of the 
low-power radio AGN out to z = 1.3. We combine our 
results with the new findings and ideas on the cosmolog- 
ical relevance of radio AGN in order to study the impact 
of radio luminous AGN on galaxy formation. 

The paper is outlined as follows. In Sec. [2] we define 
the VLA-COSMOS AGN sample. In Sec. 1 we derive the 
radio luminosity function for the VLA-COSMOS AGN, 
and extend it to hig h radio powers using the results from 
IWillott et~afl (|200lh . In Sec. Hwe constrain the evolu- 
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tion of radio AGN out to z = 1.3. In Sec. [5] we analyze 
the properties of local and intermediate redshift weak 
and powerful radio AGN; in Sec. [5] we compute the radio 
AGN mass function, and derive and compare the star for- 
mation quenching and radio- AGN triggering rates. We 
discuss our results in the context of galaxy formation and 
evolution in Sec. and summarize them in Sec. [5] 

We report magnitudes in the AB system, adopt Hq = 
70, Am = 0.3, fiA = 0.7, and define the radio syn- 
chrotron spectrum as F v oc v~ a , assuming a = 0.7. 

2. THE VLA-COSMOS RADIO AGN SAMPLE 

The sample of AGN galaxies used here is presented 
in ISmolcic et al.l ((2008a, S08a hereafter), and we briefly 
summarize it below. 

Using radio and optical data for the COSMOS field, 
S08a have constructed a sample of 601 AGN galaxies 
with z < 1.3 f rom the entire VLA-C OSMOS Large 
Project catalog ( Schinnerer et al.ll2007l ). The selection 
required optical counterparts with iab < 26, accurate 
photometry, and a S/N > 5 (i.e. > 50 /iJy) detection 
at 20 cm, and is bas ed on a rest-frame optical color 
classification (see also ISmolcic et al.l I2006D . The clas- 
sification method was well-calibrated using a large lo- 
cal sample of galaxies (~ 7,000 SDSS "main" spec- 
troscopic galaxy sample, NVSSand IRAS surveys) rep- 
resentative of the VLA-COSMOS population. It was 
shown that the method agrees well with other indepen- 
dent classification schemes based on mid-infrared colors 
(jLacv et alj |200j; [Stern et all l200j) and optical spec- 
trosc opy (|Baldwin. Phillips k Terlevichj|198lHBest et all 
2005). The selected sample of AGN is estimated to be 
~ 90% complete. 

The rest-frame color classification procedure efficiently 
selects mostly type 2 AGN (such as LINER, Seyferts), 
and absorption-line AGN (with no emission lines in the 
optical spectrum), while type 1 AGN (i.e. quasars, < 20% 
of the total AGN sample) are not included in the current 
sample (see S08a for detailed definitions of the samples). 

Although based on a color identification (as opposed 
to an optical spectroscopic classification), the output 
of the selection of our intermediate redshift (z < 1.3) 
AGN sample is comparable to those of numerous local 
(z < 0.3) radio AGN samples extensively stud i ed in 
the literature fe e. ISadler et al.l l2002t IBest et all 120051 : 
iMauch fc Sadler! I2007L drawn from the SDSS, 2dF, and 
6dF optical surveys combined with the NVSS and FIRST 
20 cm radio surveys). This is an important feature as it 
enables a straight-forward and fair comparison of our re- 
sults with those based on these local studies (e.g. the 
local AGN radio luminosity function). 

Out of the 601 selected AGN galaxies 262 have spec- 
troscopic redshifts, while the remaining sources have 

very reliable photometric redshifts available (a (^j^j = 

0.027; see S08a and references therein). Based on Monte 
Carlo simulations, S08a have shown that the photomet- 
ric errors in the rest-frame color introduce a small, ~ 5%, 
uncertainty in the number of the selected galaxies. Here 
we use the S08a sample of AGN galaxies, statistically 
corrected for this effect. 

3. THE 1.4 GHZ LUMINOSITY FUNCTION FOR RADIO 
AGN IN VLA-COSMOS 



3.1. Derivation of the luminosity function (LF) 

We derive the radio LF ($) for our AGN galaxies in 
four redshift b ins using the standard 1/Knax method 
(| Schmidt! 11968) . We limit the accessible volumes a) on 
the bright end by the minimum redshift cut-off of the 
redshift range in consideration or the minimum redshift 
out to which an object could be observed due to the 
optical saturation limit of i* = 16 (AB mag; see also 
iCapak et~a l. 2007), and b) on the faint end by the max- 
imum redshift out to which a galaxy could be observed 
given the flux limits on both the radio and optical data. 
In practice, the latter is dominated by the radio detection 
limits, rather than the optical, as the major fraction of 
the sources used here has «ab band magnitudes brighter 
than ~ 24 (see Fig. 21 in S08a). 

We further take into account the non-uniform rms 
noise level in the VLA-COSMOS mosaic via the differ- 
ential visi bility area (i.e. areal cov erage, Ak, vs. rms; see 
Fig. 13 in iSchinnerer et all 2007). Hence, for a source 
with a 1.4 GHz luminosity Lj its maximum volume is 

V max (Lj) = X}fc=l Ak ■ Knax(^max! 

In order to robustly derive the LF we take several ad- 
ditional corrections into account: a) the VL A-COSMOS 
detection completeness ([Bondi et al.l 120081) . and b) the 
AGN galaxy selection bias due to the rest-frame color 
uncertainties. We correct for these in the same way 
as described in ISmolcic et alj (|2008bl S08b hereafter). 
The median of the first correction (as a function of flux 
density) is ~ 10% (reaching a maximum value of 60% 
in only one of the l owest flux density bins; see Tab. 1 
in lBondi et al.l 12001t h The second correction, based on 
Monte Carlo simulations, yields a reduction of ~ 5% of 
the average number of radio AGN (see below; see also 
S08b for a more detailed description). 

Hence, in the i th luminosity bin the comoving space 
density (<&;), and its corresponding error (<Xj), are 
computed by weighting the contribution of each (j th ) 
galaxy by the co mpleteness correction factor, f<j e t (see 
iBondi et al.l[200l : 
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The selection bias due to the rest-frame color uncer- 
tainties is accounted for via Monte Carlo simulations. As 
described in S08a and S08b, in each iteration the rest- 
frame color error distribution is generated (see Fig. 5 in 
S08a) and added to the rest-frame color derived by SED 
fitting. AGN galaxies are then re-selected, and the LF 
is derived as described above. In this way we obtain 100 
realizations of ($^,(7^) for each luminosity bin, and we 
take the median values as representative. 

3.2. The radio AGN luminosity function 

The 1.4 GHz radio LFs for our AGN galaxies for 
the 4 chosen redshift bins are shown in Fig. [TJ and 
tabulated in Tab. [T] In each panel in Fig. [1] we 
also show the local 2 cm LFs fo r AGN derived by 
ICondon et ail (I2002D. ISadler et al l ([20021) . IBest et all 
(I2005h . and IMauch fc Sadler! (|2007ft . Our derived LF in 
the local redshift bin (top left panel) agrees well with 
the local LFs. This is quite remarkable as a) our identi- 
fication method is based only on photometry contrary to 
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TABLE 1 

Luminosity functions for VLA-COSMOS 
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the spectroscopic selection performed by the other stud- 
ies (ISadler et al.ll2002t iBest et al.ll2005t iMauch fc Sadlerl 
|2007j), and b)~the 2D° COSMOS field samples a signif- 
icantly smaller comoving volume at these low redshifts 
compared to the almost all-sky surveys used by the other 
studies. This verifies that both our selection method, 
as well as the derivation of the LF are correct. Within 
the errors our lowest-redshift LF se ems to be be s t rep- 
resented by the local LF derived by I Sadler et al.l (2002, 
Sad02 hereafter), and we adopt this local LF for further 
analysis of the evolution of our AGN population. 

In the top right panel in Fig. [T] we compare our de- 
rived AGN LF ( 35 < z < 6) wit h the volume densi- 
ties computed by ISadler et al.l ((2007) for a redshift range 
0.4 < z < 0.7 based on 2SLAQ Luminous Red Galaxy 
Survey and NVSS, FIRST data. We convert their data 
to a base of dlogL, and scale them to match the mean 
redshift of our redshift bin (z = 0.475 compared to their 
z = 0.55) using pur e luminosity evolution as obtained by 
ISadler et al.l (|2007l note that this only decreases their ra- 
dio luminosities by 10%). The LFs are in excellent agree- 
ment. It is noteworthy that we well constrain the volume 
densities at the low luminosity end (< 10 25 W Hz -1 ), 
while the sample of ISadler et al.l (|2007h extends further 
out to 1.6 x 10 27 W Hz" 1 (see also Fig. El and SecEJ). 

The VLA-COSMOS AGN volume densities in the 
two highest redshift bins (0.6 < z < 0.9 and 0.9 < 
z < 1.3) are shown in the bottom panels in Fig. [TJ 
The radio AGN LF at these redshifts has been stud- 
ied in full detail only for higher power AGN (Li.4GHz> 



10 25 W Hz" 1 : iDunlop fc Peacock! H9901 iWillott et all 
I200H iWaddimrton et all 1200 ID ." while constraining the 
low-luminosity end (sampled here) has been hampered 
by low number statistics due to small field sizes. 
ICowie et al.l ((2004) combined two deep 1.4 GHz radio 
surveys of the HDF-N (40' in diameter; 5a ~ 40 fiJy in 
the central region) and SSA13 (34' in diameter; 5cr ~ 
25 fiJy in the center) fields, and derived the AGN radio 
LF in two separate redshift ranges similar to ours, and at 
comparable luminosities to our work. Our r esults are in 
qualit ative agreement with those derived bv lCowie et al.l 
(|2004l c.f. their Fig. 3). However, their sample is a factor 
of ~ 7 smaller than the one used here. 

4. THE EVOLUTION OF RADIO AGN 

In this section we constrain the evolution of our low- 
power radio AGN using the VLA-COSMOS AGN data 
(Sec. 14. ip . We further extend this to hig h-power ra- 
dio A GN using the model obtained by IWillott et al.l 
(2001) based on high -power AGN samples drawn from 
the 3CRR, 6CE, and 7CRS surveys (Sec.EOJ. 

4.1. The evolution of low-power radio AGN galaxies in 
the COSMOS field 

The evolution of an astrophysical population is usu- 
ally parameterized by monotonic density and luminosity 
evolution of its local luminosity function: 



$ 2 (i) = (l + z) QD x $ 2 



L 



(l + z) a ^ 



(2) 



where ur> and are the characteristic density and 
luminosity evolution parameters, respectively, L is the 
luminosity, $> Z (L) is the luminosity function at redshift 
z, and $ z= o is the local luminosity function. 

It is well known that strong degeneracy between lumi- 
nosity and density evolution exists, even if the observa- 
tional data sample the turn-over (' knee') of the LF at dif - 
ferent cosmological times (see e.g. lLe Floc'h et aL¥ 2005). 
The VLA-COSMOS AGN sample in particular lacks high 
luminosity objects that could constrain the turn-over of 
the LF in all our redshift bins (see Fig. [T|). Hence, we 
do not attempt to try to break the luminosity - density 
evolution degeneracy, but we separately constrain both 
pure density (PDE; ar, = 0) and pure luminosity (PLE; 
old = 0) evolution based on our data. We adopt the 
Sad02 local AGN LF ($2=0) as the representative LF in 
the local universe, which is given by the following ana- 
lytical form: 



" L ' 


l— a 






exp i 


2a 2 



log(l + A) 



(3) 



where a = 1.58, a = 1.0, $* = 7.6 x 10~ 6 Mpc~ 3 , 
and L* = 2.1 x 10 24 W Hz" 1 for their AGN population 
(scaled to the cosmology used here, and to the base of 
dlogi). 

For both PDE and PLE we derive the evolution by 
summing the x 2 distributions obtained for a large range 
of evolution parameters and [(a^fi) for PLE; 
(0,ar>) for PDE] in each particular redshift bin (exclud- 
ing our first - local - redshift bin). The uncertainty in 
a is then taken to be the la error obtained from the 
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Fig. 1. — 1.4 GHz luminosity functions (LFs) for AGN in the VLA- COSMOS survey (fil led red circles), show n for four redshift r anges 
indicated in each panel . The local AGN volume densities derived by [Co ndon ct al. (2002], dashed black line), ISadler et al.l 1120021 , solid 
orange line), [Best et al] (J2005J, empty squares), and[Mauch Sz Sadler (2007, dotted line) are also plotted in each pane l . The empty symbols 
in the top right panel are the radio AGN volume densities derived for the redshift range 0.4 < z < 0.7 by [Sadler ct al. (2007J) for the 2SLAQ 
Luminous Red Galaxy survey sample, and scaled here to match our mean redshift value of 0.475 (see text for details). 



X 2 statistics. Our results yield a pure density evolution 
with old = 1.1 ±0.1, or alternatively a pure luminosity 
evolution with a_L = 0.8 ± 0.1. 

In Fig. [5] we show the luminosity density for our AGN 
in the four redshift ranges defined in Fig. [TJ We also 
plot the luminosity density given the derived PD and PL 
evolutions (lines). From the figure it becomes obvious 
that we cannot distinguish between these two types of 
evolution given our data, as our sample, tracing only the 
lower luminosity end of the luminosity function, appears 
to be described equally well by both. Thus, in the fur- 
ther analysis we will take the range of the two best fit 
evolution models (taking also their errors into account) 
as representative of the range of uncertainties. 

In summary, our results imply either a pure luminos- 



ity evolution with L» oc (1 + ^) u - 8±01 or a pure density 
evolution where <&» oc (1 4- z) 11±01 . Regardless of which 
model is physically more appropriate to describe the real 
cosmic evolution of the VLA-COSMOS AGN population, 
both imply a modest evolution of radio luminous AGN 
in the luminosity range of ~ io 22 ~ 25 W Hz -1 . 

Our results are consistent with previous findings. 
B ased on a V/Vmax analy sis out to a redshift of ~ 
1, IClewlev fc Jarvia ()2004f ) have found no strong evo- 
lution of low radio luminosity sources (L325MHZ ^ 
10 25 W Hz -1 sr _1 , corresponding to Li.4ghz J$ 4.5 x 
10 25 W Hz" 1 ) a t least up to z ~ 0.6. Further, 
I Sadler et al.1 (|2007f ) have constrained the evolution of the 
luminosity function for radio luminous AGN (Li 4GHz ~ 
10 24-27 w Hz -i) out to z _ 7 Tj sing t i ie i oca i AGN 
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Fig. 2.— Luminosity density for the VLA-COSMOS AGN in four 
redshift ranges (filled circles). The lines represent the best fit pure 
density (dashed line) and pure luminosity (solid line) evolution of 
the local 20 cm AGN luminosity function (see text for details). 

LF given by iMauch fc Sadler! (|2007l ) they have found 
a pure luminosity evolution with oc (1 + 2;) 2 0±0 - 3 . 
The VLA-COSMOS AG N LF in th e reds hift range 
0.35 < z < 0.6, and the ISadler et all p007h AGN LF 
(0.4 < z < 0.7) agree remarkably well (see top right 
panel in Fig. Q] & Fig. [3]) . If we parameterize the evo- 
lution of the VLA-COSMOS AG N in the same way as 
described in ISadler et al.l (|2007l using the same local 
LF), we derive a pure luminosity evolution of 1.2 ± 0.4 
for the 0.35 < z < 0.6 redshift bin. It has to be 
noted, however, that the VLA-COSMOS and the 2SLAQ 
samples constrain different radio luminosity ranges (see 
Fig. [3]). Our results agree withi n the uncertainties with 
those of IMauch fc Sadler] (|2007l ), however they on aver- 
age imply a slower evolution. This is consistent with 
a more modest evolution of low ra dio-power, compared 
to high radio-power, radio AGN (IWillott et all 120011 : 
IWaddington"etai1 l200lt IClewlev fc Jarvisl I2004D . Our 
findings are also in very good agreement with the recent 
results based on NVSS/FLRST and the Me gaZ-Luminous 
Red G alaxy catalog drawn from the SDSS (jDonoso et alJ 
l2008h . 

4.2. The evolution of high-power radio AGN 

In order to quantify the contribution of different popu- 
lations to the comoving radio energy density at different 
cosmic times, all radio AGN populations need to be con- 
sidered. As low- and high- power radio AGN seem to 
evolve in a different manner (see Sec. [IJ and our VLA- 
COSMOS AGN sample only the low power radio AGN, 
we need to make the following assumptions about the 
evolution of the high power radio sources. 

Based on the 3CRR, 6CE, and 7CRS radio sur- 
veys combin e d wi th complete optical spectroscopy, 
IWillott et alJ (|2001h have successfully modeled the radio 



AGN LF using two radio populations - a less powerful 
(Li.4GHz ~ 2.5 x io 25 ~ 27 W Hz -1 ) population compris- 
ing both FR I and FR II sources, and a powerful popu- 
lation (Li.4GHz > 2.5 x 10 26 W Hz -1 ) comprising mostly 
FR II sources. They have modeled the evolution of the 
first population as a pure density evolution up to a maxi- 
mum redshift (~ 0.7) beyond which any evolution ceases. 
The evolution of the powerful population has been as- 
sumed to change in density following a Gaussian distri- 
bution in redshift, which was allowed to have a different 
shape beyond its r edshift peak at z ~ 2 (see Tab. 1 in 
IWillott et a!1l200lD . It is wort h noting that the Wi l lott et 
al. model agrees well with the lDunlop fc Peacock! (|1990l ) 
steep spectrum model. 

The VLA-COSMOS AGN sample constrains the faint 
end of t he radio AGN lumin osity function, and here we 
use the IWillott et alJ ( |200lf ) model to extend our radio 
LF to high powers. In Fig. [3] we compare our LFs with 
the IWillott et a l. (2001) model, after the latter has been 
converted to the current cosmology and the 151 MHz ra- 
dio luminosities scaled to 1.4 GHz. The VLA-COSMOS 
AGN data and the Sad02 LF constrain the radio LF more 
robust ly at the faint end compared to the IWillott et al.l 
(|2001f ) model for their less-powerful radio AGN. Thus, in 
the further analysis we will constrain the low-power ra- 
dio AGN LF and its evolution using the VLA-COSMOS 
sample as described in the previous sectio n, and we will 
use on ly the powerful radio AGN model by IWillott et al.l 
(|2001h to describe the evolution of high radio-power AGN 
(see dash-dotted curves in Fig. [3]). 

4.3. The evolution of the comoving radio luminosity 
density for A GN galaxies 

At a specific cosmic time the integrated comoving lumi- 
nosity density represents the total power per unit comov- 
ing volume of a given astronomical population. Thus, 
if divided into distinct populations of objects it traces 
their relative contribution to the overall power output 
at a given redshift. To estimate the contribution of low 
and high radio power AGN to the AGN radio luminos- 
ity output at different cosmic times, we investigate in 
the following the evolution of the comoving 1.4 GHz lu- 
minosity density, 4G Hz, for our VLA-COSMOS AGN 
(Li. 4G Hz < 10 25 W Hz^ 1 ) and the hig h-power AGN 
(kuLGHz > 10 26 W Hz" 1 ), adopting the IWillott et all 
(2001) model. 

For a given redshift ^i.4GHz has been computed by 
integrating the luminosity density over the entire range 
of radio luminosities. For the VLA-COSMOS AGN the 
luminosity density has been constrained using our best 
fit PLE and PDE models (see curves in Fig. [2] and 
Sec. 14. 1|). and for t he high-power population using the 
IWillott et alJ (|2001f) model 'C, with t he corresponding 
errors (see tab. 1 in lWillott et~aTll2001l and Sec. 021). 

In Fig. |4] we show the evolution of the comoving 20 cm 
integrated luminosity density for all radio AGN as well 
as the low and high radio power AGN separately. The 
evolution of these two populations is very different; high- 
power AGN evolve significantly stronger than low-power 
AGN. The total fii. 4 GHz for radio AGN is dominated by 
low-power AGN at low redshifts (z < 0.7) where the 
contribution of high-power radio sources is negligible. 
However, at z > 0.7 high-power sources begin to con- 
tribute significantly to the overall integrated luminosity 
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Fig. 3. — VLA-COSMOS AGN volume densities at 20 cm in four redshift ranges (red filled circles; analogous to Fig.flJ. Also shown in 
each panel is our best fit evolution corresponding to the range given by pure density and pure luminosity evolution o f the Sad02 local LF 
and their corresponding errors (orange shaded curve), as well as the radio AGN luminosity function model given by IWillott et alj J200ll , 
model 'C') for their less luminous population (dashed blue curves), and high-luminosity population (dash-dotted blue curv es; scaled to 
curren t cosmology and 1.4 GHz radio frequency; see text for details). In the top right panel the volume densities derived bv lSadler et al.l 
(120071 . scaled to match our mean redshift value of 0.475; see also Fig. [TJ arc shown. 



density, and at a redshift of ~ 1.3 their contribution to 
the total AGN £li.4GHz is comparable to the one of low- 
power AGN. The implications of these populations for 
galaxy formation and evolution, also out to higher red- 
shifts (z = 2.5), are discussed in Sec. [7] 

5. PROPERTIES OF RADIO AGN 

In this section we outline and compare the properties 
of low and high radio-power AGN in the local (Sec. 15. ip 
and intermediate-redshift (Sec. I5.2|) universe. We find 
that already by z ~ 1 the host galaxies of low-power, 
VLA-COSMOS AGN have built-up stellar and black hole 
masses comparable to the highest mass galaxies observed 
locally. As their black hole masses are already signif- 
icant at these intermediate redshifts, this implies that 



they will no longer be able to have a phase of high accre- 
tion (i.e. vigorous BH growth), consiste nt with numerous 
previous studies of similar samples (e.g.lAllen et al.ll2006t 
lEvans et al]l2006t lHardcastle et al1l2006Ll2007t ). For the 
high radio power AGN on the other hand, extensive ev- 
idence exists in the literature that they accrete at high 
rates - representing a mode of substantial BH growth. 

5.1. Local universe 

Various correlations are found in the literature between 
the presence of emission lines in AGN and their e.g. ra- 
dio power, black hole and stellar mass, as well as envi- 
ronment and the galaxies' gas content. We outline these 
below. 

First, almost all FR I - low power - radio galaxies are 



V. Smolcic et al. 




z 



Fig. 4. — Evolution of the comoving 20 cm integrated luminos- 
ity density for VLA-COSMOS AGN (orange curve) galaxies since 
z = 1.3. Shown is also the evolution of th e high-luminosity ra- 
dio AGN, adopted from IWillott et all ||2001| . hatched region; the 
thick and dashed lines correspond to the mean, maximum and 
minimum results, respectively). The evolution for the total AGN 
population, obtained by co-adding the VLA-COSMOS and high- 
luminosity AGN energy densities, is shown as the red-shaded curve 
(see text for details). 

LERAGN, while optical hosts of FR lis, which are typi- 
cally more p owerful than FR Is jf k naroff fc Rileyl ll974b 
IOwenlll993 : Lcdlo wfe OwerJfl996| ). usually have strong 
emission lines 18 . Recently, based on a large statistically 
significant sam ple of local radio - optica l sources (SDSS- 
NVSS-FIRSTV lKauffmann et all ((20081) have found that 
the fraction of radio AGN whose optical hosts have emis- 
sion lines in their spectra (predominantly HERAGN) is 
a strong function of radio luminosity. This emission- 
line galaxy fraction is roughly constant (~ 40%) up to 
Li.4GHz ~ 10 25 W Hz -1 , beyond which it steeply rises 
a pproaching ~ 80% at ~ 4 x 10 25 W Hz" 1 (see Fig. 3 
in Kauffm ann et al.ll2008f ). This 'critical' luminosity, ob- 
served by Kauffmann et al., is remarkably close to the 
FR I - II break luminosity, as well as to the power 
which roughly separates the radio sources which show 
strong cosmological evolution from those which do not 
(see e.g. Fig. 2] and Seed]). Furthermore, based on the 
results of Kauffmann et al., the luminosity o/Li.4qhz ~ 
10 25 W Hz^ 1 can be thought of as a rough threshold be- 
tween high- and low- excitation radio AGN. Thus, most 
of the low radio power AGN (Li. 4G Hz < 10 25 W Hz" 1 ) 
are LERAGN, while the majority of powerful radio AGN 
(Li.4GHz > 10 25 W Hz" 1 ) are HERAGN. 

On the other hand, the fraction of emission line (i.e. 
high-excitation) radio AGN in the local universe strongly 
decreases as a function of both stellar mass and v elocity 
dispersion (see Fig. 3 in iKauffmann et al.l l2008t ). This 
implies that, at least locally, high-excitation (or alterna- 
tively high radio power) AGN tend to have lower stellar 

18 Note, however, that the correspondence between the FR class 
and the presence of emission lines is not exactly one-to-one. 



masses, as well as lower black hole masses compared to 
LERAGN. Even further, the latter constitute the most 
massive galaxies in the universe (Mq > 10 11 M Q ) that 
preferentially occupy the centers of high galaxy density 
regions (Baum et al. 1992, Best et al. 2005). 

Furthermore, various studies in the literature have 
shown that HERAGN, contrary to LERAGN, tend 
to show unusually blue off-nuclear continuum colors 
( Baum et al.1 Il992f). ev idence of recent star formation 
(Baldi fc Capettir [2008) . and often have distorted opti- 
cal morphologies sugg esting that they have undergone a 
recent major merge r (Hcck man et al.lfl986t [Baum et al.1 
119921: iBaldi fc Capettill2008ft . CO and HI observations 
of radio galaxies suggest larger amounts of cold gas in 
powerful, compared to low-power, radio sources (see e.g. 
Fig. 8 in Evans et al. 2005, see also Emonts et al. 2008), 
and there is evidence supporting qu antitatively larger 
amounts of dust, and therefo re gas (|Leon et al.l 120011 : 
ISolomon fc Vanden Bout] l2005f ) in HERAGN compared 
to LERAGN (de Kofi et al. 2000, Miiller et al. 2004). A 
summary of the properties of radio AGN in the nearby 
universe is given in Tab. [2j 

5.2. Intermediate redshift 

In this section we investigate the properties of low and 
high radio-power AGN at intermediate redshifts, and 
compare these to the properties of their local counter- 
parts. In Sec. [3] and Sec. 2] we have derived the 1.4 GHz 
radio luminosity function, and its evolution, for a sam- 
ple of ~ 600 radio luminous AGN (0.1 < z < 1.3) 
drawn from the VLA-COSMOS survey. One of the 
main characteristics of this radio AGN sample is that 
it consists of predominantl y (96%) low-powe r AGN with 
Li.4GHz < 10 25 W Hz" 1 . ILedlow fc Owen] (fl996l) have 
shown that the threshold between weak FR I types of 
galaxies and powerful FR lis is a function of the host 
galaxy's optical (R band) magnitude, Mr. Hence, to bet- 
ter constrain the properties of the VLA-COSMOS AGN 
in Fig. [5] we plot them in the Li.4GHz vs. Mr plane. Al- 
most all (98%) of our AGN occupy the weak radio AGN, 
FR I region in this plane. 

In order to get a further insight into the composition 
of the VLA-COSMOS AGN, we ma ke use of the sources 
with available optical spectroscopy. iTrump et alj (|2007l ) 
have carried out a spectroscopic survey of X-ray and ra- 
dio selected AGN candidates in the COSMOS field (using 
the Magellan/IMACS instrument). They have classified 
their sources into a) broad emission line AGN ('bl'), b) 
narrow emission line AGN ('nl'), c) red galaxies with only 
absorption features in their spectrum ('e') and d) hybrid 
objects showing a mix of narrow emission lines and a 
red galaxy continuum shape ('nle'). Combining a radio 
AGN detection with this classification scheme, we can 
very roughly consider the red and hybrid galaxy spec- 
troscopic classification ('e'+'nle') as a proxy for LER- 
AGN, and the narrow emission line AGN classification 
('nl') as a proxy for HERAGN. About 35% of our VLA- 
COSMOS AGN have an available classification based 
on IMACS spectroscopy, and in Fig. [6] the fractions of 
'e+nle' and 'nl' galaxies are shown as a function of both 
redshift (left panel) and radio luminosity (right panel). 
The absorption line and hybrid galaxies dominate the 
VLA-COSMOS sample at a constant ~ 80% level at all 
redshifts [z < 1.3), while the narrow emission line AGN 
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TABLE 2 

Properties of LERAGN and HERAGN in the local universe 





LERAGN 




HERAGN 




Property 




Ref. 




Ref. 


Object class 


mainly FR I 


(1),(2) 


mainly FR II 


(1),(2) 


radio luminosity 


Ll.4GHz < 10 25 W Hz' 1 


(3) 


Ll.4GHz > 10 25 W Hz" 1 


(3) 


density environment 


moderate-to- high 


(4),(5) 


moderate-to-low 


(4) 


optical morphology 


regular 


(8),(9) 


often distorted 


(8),(9),(10) 


optical color 


red 


(9),(H) 


bluer (compared to LERAGN) 


(10),(9) 


stellar mass 


highest (> 5 X 10 10 M Q ) 


(3),(H) 


lower (compared to LERAGN) 


(3),(12) 


BH mass 


highest (~ 10 9 M Q ) 


(9),(13) 


lower (compared to LERAGN) 


(3),(14) 


ISM content 


low 


15),(16),(17) 


higher (compared to LERAGN) 


(15),(l6),(17),(18),(19) 


accretion mode 


radiativcly inefficient 


(1),(2),(14),(20) 


radiativcly efficient 


(1),(2),(20),(21),(22) 



References: fl) lEvans et al,l 120061. 



Note. 

(5) IBest et Id] 120051 , (8) IHeckman et al.111989, (9)IBaum et al.l~992l, (lO)IBaldi fc Capettil 120081. (ll)ISmolcic et al.M2008al, (12) 



I Tasse et all 12008. f 13) this work (14) IChiaberge et a~l2(~i5l, (15) ILeon et al. 1120011 . ( 16) Ide Koff et al.l 11999. (17) IMiiller et al.l 
2004, (18) Evans et al. 2Q0~, (19) lEmonts et~l . 2008, (20) Hardcastlc e t al.M2007l. ~f21) Barthel 1989, (22) Haas 2004 
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Fig. 5. — Monochromatic 1.4 GHz radio power for VLA- 
COSMOS AGN as a function of their host galaxy absolute R band 
magnitude. The dashed line corresponds t o the separation betwee n 
FR I and FR II types of galaxies given by|Ledlow & Owen (1996). 
Note that almost all of the VLA-COSMOS AGN occupy the low- 
power FR I region of this plane. 

contribute about ~ 20% at all redshifts. There is an in- 
dication that in the highest redshift bin (1 < z < 1.3), 
where our most luminous (Li.4GHz > 10 25 W Hz- 1 ) ra- 
dio AGN are observed (see Fig.QJ, the fraction of narrow 
emission line AGN rises. However, given the large error- 
bars in this redshift range (due to the low number of 
spectroscopically observed radio AGN) we cannot draw 
any robust conclusions regarding this. In the right panel 
of Fig. [6jwe show the fraction of 'e+nle' and 'nP galaxies 
as a function of 1.4 GHz power. Although the number of 
sources with Li.4GHz > 10 25 W Hz -1 is low in the VLA- 
COSMOS AGN sample, there is an indication that the 
fraction of narrow-line objects (thus roughly HERAGN) 
increases beyond ~ 10 25 W Hz" 1 , consistent with the 
findings in the local universe (|Kauffmann et al.ll200~h . 

Assuming that the spectroscopic sub-sample represents 
well the full sample (note that this is a rather robust as- 
sumption as shown in S08a; see their Figs. 3 and 21), we 
conclude that LERAGN dominate the VLA-COSMOS 



radio AGN sample. Thus, similar to the findings in the 
local universe, our low radio power AGN at intermediate 
redshift are preferentially LERAGN . 

Stellar masses, using a Chabrier (2003) initial mass 
function (IMF), have been computed for the entire VLA- 
COSMOS galaxy sample (z < 1.3) by S08a. In Fig. [7] 
(top panels) we show the stellar masses for our VLA- 
COSMOS AGN as a function of redshift and 1.4 GHz 
luminosity. The median stellar mass of our radio AGN 
is 1.6 x 10 11 M , i.e. < logAf* >= 11.2 (in any given 
redshift range) with a la scatter of 0.4 dex. This is 
consistent with th e stellar masses of the most massive , 
local galaxies (e.g. iBaldrv et all 12004 IBest et alJ[200~J) . 
In addition, it is worth noting that their average rest- 
frame colors are consistent with red galaxy colors (see 
e.g. Fig. 9 in S08a). 

We further compute the BH masses for the full VLA- 
C OSMOS AGN sample us ing the local correlation given 
bv lMarconi fc Hurra (|2003f ) which relates the K-band lu- 
minosity to the BH mass: 

log 10 M BH = 8.21 + 1.13 x (log 10 L K - 10.9) (4) 

Lk the rest-frame K-band luminosity (also in so- 
lar units). The scatt er in the relation is 0.5 dex 
(|Marconi fc HuntJ [2003). The K-band rest-frame lu- 
minosity for our VLA-COSMOS AGN was computed 
via SED fitting as described in detail in S08a. In or- 
der to take into a c count passive luminosity evolution 
(see iHopkins et al.l l2006bl and references therein) we 
change the normalization constant of the above rela- 
tion as a func t ion of r edshift following the results from 
Hopk ms et~aLl (|2006U see short-dashed line in the bot- 
tom left panel in their Fig. 2). In Fig. [7] (bottom panels) 
we plot the estimated BH masses as a function of both 
redshift and radio 1.4 GHz luminosity. The median BH 
mass of our radio AGN is log 10 Mbh = 8.8 with a stan- 
dard deviation of 0.4. Note that the apparent trend of 
BH mass with radio luminosity (bottom right panel in 
Fig. [7J) reflects the dependence of radio luminosity on 
redshift in our radio flux limited sample (see Fig. 16 in 
S08a) rather than a real trend. 

The black hole masses of our intermediate redshift ra- 
dio AGN are comparable to the highes t black hole masses 
known in t he local universe (see e.g. iMcLure fc Dunlopl 
120021 12004f ). This implies that the low -power radio AGN 
have assembled their BH masses already by these in- 
termediate redshifts, and that their BHs cannot grow 
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Fig. 6.— The fraction of VLA-COSMOS AGN with available Magellan/IMACS spectroscopy IITrump et alj|2007h as a function of 
redshift (left panel) and 1.4 GHz radio luminosity (right panel). The galaxies have been divided into two sub-samples, i) narrow emission 
lin e (nl) galaxies, and ii) elliptical galaxies (e) plus composite objects showing a red galaxy continuum and narrow emission lines (nle; 
sec Trump ct al. 2007 for details). These two subsamples roughly correspond to high- and low-excitation AGN, respectively (see text for 
details). The quantity and symbols for objects of each type is indicated in the left panel. Note that there is a slight (although noisy) trend 
showing that the fraction of narrow emission line objects (i.e. roughly high-excitation AGN) increases as a function of both redshift and 



radio luminosity. 
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Fig. 7. — Logarithm of stellar (top panels) and black hole 
(bottom panels) masses for the VLA-COSMOS AGN as a func- 
tion of redshift (left panels) and 1.4 GHz radio luminosity (right 
pan els). The stellar masses w ere computed via SED fitting using 
the Bruzual & Chariot (2003) stellar population synthesis models 
jChabrierl 120031 IMF; see S08a for full details). The BH masses 
were derived using eq. [4] and corrected for passive luminosity evo- 
lution (see text for details). The dashed line in the right column 
corresponds to the lowest 1.4 GHz luminosity observable out to 
z = 1.3 given the VLA-COSMOS Large Project detection limit of 
~ 50 ^Jy beam -1 . 

significantly since z ~ 1.3. Therefore, they must be in 
a mode of modest growth in BH mass. This result is 
in agreement with numerous findings in the literature 
implying that LERAGN accrete radiatively inefficiently, 



at sub-Eddington rates (I Allen et all l2006t I Evans et al.l 
2006; lHardcastle et"aT1l2006L 120071 

Given all of the above, the composition of our VLA- 
COSMOS AGN, the majority of which are shown to be 
LERAGN, is consistent with the properties of low radio- 
power AGN in the local universe: Our intermediate red- 
shift AGN have already by z ~ 1 assembled both their 
stellar and black hole masses, comparable to the high- 
mass end of the galaxies known today. 

On the other hand, powerful radio galaxies at high red- 
shifts tend to show strong emission lines in their optical 
spectra (Rawlings et al. 1989, Baum & Heckman 1989, 
Rawlings & Saunders 1991; Willott et al. 1999, 2000), 
and are often associated with ongoing star-formation 
(Archibald et al. 2001, Greve et al. 2006, Seymour et al. 
2008) as suggested by their bluer rest-frame colors (com- 
pared to red-sequence galaxies). Again, this is consistent 
with local observations of su ch sources. In numerous 
studi es in the literature (e.g. lBarthellll989t lHaas et all 
l2004f ) these objects have been shown to accrete radia- 
tively efficiently at high (~ Eddington) accretion rates. 
Thus these HERAGN present a mode of significant BH 
growth - unlike the low radio power LERAGN. 

To summarize, the properties of local and 
intermediate-redshift radio AGN, as shown above, 
suggest that radio activity is triggered in similar popula- 
tions of objects independent of their redshift, i.e. cosmic 
time. This is in agreement with the finding that the 
rest-frame optical colors of r adio source host g alaxies 
do not change wit h redshift (jBarger et al.l 120071 S08a, 
iHuvnh et aTfl2008fl . In addition, there is converging 
evidence that low radio power AGN reflect a modest - 
radiatively inefficient - mode of BH growth, while high 
radio power AGN are undergoing a phase of significant 
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radiatively efficient - BH growth. 



6. 



THE STAR FORMATION QUENCHING AND 
RADIO-AGN TRIGGERING RATES 



One of the aims of this work is to study the link be- 
tween radio AGN activity and star formation from an 
observational perspective. Thus, in order to investigate 
whether the processes responsible for the build up of the 
red galaxy sequence and the triggering of radio AGN 
mode are related or not, in this section we study the 
stellar mass properties of the red parent (Sec. 16. lj) and 
radio AGN samples (Sec. 16. 2j) . and derive and compare 
the rates for the star formation quenching (Sec. 16. 3[) and 
(low radio power) AGN triggering (Sec. I6.4[) . We find 
no evidence that these two processes are related. This 
however does not exclude the possibility that the radio 
AGN phase may be responsible for preventing already 
assembled massive galaxies to grow higher in mass, as 
explored in Sec. 17.31 

6.1. The parent galaxy sample 

In order to derive the rates for star formation quench- 
ing and radio- AGN triggering we need to define a control- 
parent galaxy sample of our VLA-COSMOS radio AGN. 
We define this sam ple using the COS MOS photomet- 
ric redshift catalog (jllbert et alj 120081 ) from which we 
select galaxies with the same optical magnitude, red- 
shift and color criteria applied to our VLA-COSMOS 
AGN (see Sec.[2])- The rest-frame optical color PI (con- 
strained from the 3200 — 5800 A range), which was used 
to select our VLA-COSMOS AGN is not available for 
the galaxies in the COSMOS photometric redshift cata- 
log. Nonetheless, as the NUV-NIR galaxy S EP has been 
shown to be a one-param eter family (e.g. lObric et al.l 
[20061: ISmolcic et al.lf2006h . we can safely utilize another 
rest-frame color for the selection of the parent sample. In 
Fig.[8]we show the correlation between U—B and PI rest- 
frame colors for all VLA-COSMOS galaxies {z < 1.3), 
and derive an analytic relation correlating PI and U — B 
(indicated in Fig. [8|). The criterion of PI > 0.15, used 
to select the VLA-COSMOS AGN sample, corresponds 
to U — B > 0.71. Thus we select the red galaxy parent 
sample by requiring I auto < 24, 0.1 < z phot < 1.3, and 
U - B > 0.71. The U - B distribution for all (red and 
blue) galaxies, that satisfy the magnitude and redshift 
criteria, is shown in Fig. [9l Note that the adopted U — B 
threshold, based on the comparison with the PI values, 
corresponds almost exactly to the U — B value which 
separates the red from the blue galaxies. Our selection 
yields 21,525 galaxies in the red galaxy control-parent 
sample. 

6.2. The stellar mass properties of VLA-COSMOS 

AGN 

Using the 1/V max method we derive the stellar mass 
function (SMF) in two redshift bins (0.2 < z < 0.7 
and 0.7 < z < 1.3) for both our VLA-COSMOS AGN 
(Li.4GHz^ 4 • 10 23 W Hz -1 which roughly corresponds 
to a radio volume limited sample), and the red parent 
galaxy sample. The stellar mass functions are shown 
in Fig. [TUJ The SMF for the red galaxy sample agrees 
well with the re cent results based on the GOODS-MUSIC 
galaxy sample ([Fontana et al.l 12006). Similarly to what 
has been shown in the local universe for radio luminous 
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drawn from the COSMOS 
photometric redshift catalog. The red galaxy parent sample of our 
VLA-COSMOS radio AGN is selected by applying U - B > 0.71 
(dashed line) to this distribution (this is equivalent to PI > 0.15 
that was used to select our radio AGN; see Fig. 18) , 



AGN (|Best et al.l [2005), at intermediate redshifts the 
SMFs of our radio AGN sample are strongly biased to- 
ward high stellar masses. Interestingly, the redshift evo- 
lution of the SMFs of the parent galaxy and radio AGN 
samples is reversed. While for a given stellar mass the 
comoving number density of red galaxies decreases with 
redshift, the number density of radio luminous AGN in- 
cr eases. Our r esults are in good agreement with those 
of iTasse et al.l (|2008| ) based on radio observations of the 
XMM-LSS field. 

We compute the fraction of radio luminous AGN as 
the ratio of the above derived mass functions in a given 
stellar mass bin. In Fig.[TT]we show the radio AGN frac- 
tion as a function of stellar mass in the two adopted red- 
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4 ■ 10 23 W Hz -1 ; bold symbols) in two redshift ranges, indicated in the panel. For comparison, the stellar mass function of the red, parent 
galaxy sample in the same redshift ranges is also shown. 



shift ranges, and compare it to t he results based on the 
local SDSS and NVSS surveys (jBest et al.|[200l . The 
radio AGN fraction in our lower redshift (0.2 < z < 0.7) 
bin agrees remarkably well with the local findings imply- 
ing insignificant or absent evolution of the radio AGN 
fraction out to z ~ 0.7 at all stellar masses. However, 
there is a significant change in the fraction of radio lu- 
minous AGN at higher redshifts (0.7 < z < 1.3), in par- 
ticular for host galaxies with stellar masses lower than 
logM* ~ 11.3. We investigate this further in Sec. 16.41 



6.3. The star formation quenching rate 
Following Bundy et al. (2008), we define the star for- 
mation quenching rate, Q, as the fraction of all galaxies 
in a given stellar mass bin that migrate to the red se- 
quence per Gyr. This migration can occur through pro- 
cesses such as mass build up (via 'wet' mergers) or gas 
consumption (fading of an already massive galaxy to red 
colors). In Fig. [12] (left panel) we plot the fraction of red 



galaxies (as defined in Sec. 16. 1[) relative to the number of 
all galaxies (see Fig. [9]) as a function of cosmic time for 
various stellar mass bins. The red galaxy fractions have 
been derived in the same way as the radio AGN fraction 
(see Sec. I6.2j) . but using the red and total galaxy sam- 
ples. The the change in radio fraction per Gyr (a) of this 
red-galaxy-sequence build-up with cosmic time, shown in 
the left panel of Fig. [121 is then equivalent to the star 
formation quenching rate. We find that as stellar mass 
increases from logM* of 10.6 to 11.7 M Q decreases 
from 6 % Gyr" 1 to 0.6 % Gyr" 1 . Note also that for the 
highest stellar masses the red sequence is already almost 
fully in place by z ~ 1. 

6.4. The triggering rate of low radio-power AGN 

In Sec. 16. 31 we have derived the star formation quench- 
ing rate based on the build-up of the red galaxy popula- 
tion as a function of cosmic age (left panel in Fig. [T2"]) . 
In the right panel in Fig. [12] we now show the radio vol- 
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ume limited (log Li.4ghz > 23.6) fraction of our VLA- 
COSMOS AGN relative to the abundance of the red par- 
ent galaxies as a function of cosmic time, for the various 
stellar mass ranges as shown in Fig. [TTJ Contrary to the 
red galaxy fractions, our results yield that the radio AGN 
fraction decreases with cosmic age. At all stellar masses 
(10.6 < logM* < 11.7) a larger fraction of radio lumi- 
nous AGN is present at earlier, compared to later, cosmic 
times. Furthermore, the fraction of radio luminous AGN 
decreases more steeply with cosmic time for higher stel- 
lar masses. For example, at log A/* ~ 11.65 the rate of 
decrease is 10.8 % Gyr -1 , while at logM* ~ 10.75 it is 
only 0.1 % Gyr" 1 . 

Given the fraction of radio AGN at a given cosmic 
time, we can estimate the average time that each single 
massive galaxy spends as a radio AGN, r ra dio- Assum- 
ing that massive galaxies, that form the parent sample 
of our radio AGN, form around z ~ 2.5 — 3 and survive 
to z ~ 0, their lifetime is then ~ 10 Gyr. Multiplying 
this lifetime by the fraction of observed radio AGN then 
yields an approximate estimate of T ra dio- At a given cos- 
mic time, Tradio, which is in the range of roughly 10 Myr 
to 5 Gyr for our VLA-COSMOS AGN (see right panel in 
Fig. [T2"]) . is thus related either to the average duration of 
a single radio-episode or to the frequency of radio mode 
re-triggering. 

Independent calculations of radio source lifetimes, 
^radio, result in lifetimes of a few times 0.001 Gyr to a 
few times 0.1 Gyr (Alexander & Leahy 1987; Shabala 
et al. 2008). Assuming that the radio source lifetime 
does not significantly change with cosmic time, it is most 
likely that the computed r ra di reflects the occurrence of 
multiple radio-phases of a single massive galaxy at a cer- 
tain cosmic time. A radio 'on' phase of 0.001-0.1 Gyr, 



combined with the derived r ra di , implies then that ra- 
dio AGN activity has been triggered more often in the 
past than it is today, and that at a given cosmic time it 
is more often triggered in massive galaxies compared to 
lower mass galaxies. 

The minimum radio AGN triggering rate can be ob- 
tained by assuming the maximum i r adio of 0.1 Gyr, and 
taking into account the above derived r ra dio (0.010 — 
5 Gyr). It is then 0.1 % Gyr" 1 to 50 % Gyr" 1 for 
the lowest and highest stellar mass bins, respectively. 
Yet, in the previous section we have computed a star 
formation quenching rate of 6 % Gyr -1 for the lowest 
and 0.6 % Gyr -1 for the highest stellar mass bins (see 
left panel in Fig. [T2")) . Hence, comparing these estimates 
yields that our derived radio AGN triggering rates are 
by many factors off the derived star formation quench- 
ing rates, especially for the highest stellar masses. 

If the quenching of star formation, that causes the 
build-up of the red galaxy sequence, is related to a single 
episode (i.e. once the star formation is suppressed it does 
not restart at a significant level), then it is highly unlikely 
that a direct connection between the two phenomena - 
the low radio power AGN triggering and the star for- 
mation quenching rates - exists. However, these results 
do not exclude the possibility that low radio power AGN 
are responsible for preventing the galaxies to grow higher 
in mass, once they have been established as the most 
massive 'red-and-dead' galaxies, as has been proposed 
in numerous cosmological models (ICroton et all 120061: 
i Bower et afll2006t ISijacki fe SpringeW2006t iSijacki et all 
[2007D- This will be further explored in Sec. El 

7. IMPLICATIONS FOR GALAXY FORMATION 

7.1. High and low radio power AGN: Different stages 
of galaxy evolution 

In the last decades studies at various wavelength 
regimes have converged towards a widely accepted galaxy 
formation picture. Galaxies are thought to evolve 
in time from an initial stage with spiral morphology 
and blue optical colors towards elliptica l morphologies 
with red optical c olors (iBell et al.ll2004a|lbl: iBorch et al] 
2006tlFaber et al.l2007HBrown et al.l2007tlHopkins et all 



120071 ). This evolution is not linear but happens through 
interspersed episodes of massive mass accretion onto 
the stellar body as w ell as their central mass i ve black 
holes seen as quasars ((Sanders fc Mir abel 1996; Sanders 
2003). The blue-to-red galaxy evolution is accompa- 
nied with significant build-up of stellar and BH mass; 
the reddest galaxies observed in the universe are also 
the oldest and the mo s t massive ones (My > 1 11 M Q ; 
e.g. iBaldrv et all 12004 [20061 : iFaber et al.ll2007h . Merg- 
ers between galaxies, rich in cold gas ('wet' merg- 
ers), appear to be one of the key processes governing 
this mas s build-up and the blue-to-red galaxy transfor- 
mation (ISanderd [20031 IBell et all l2004allbl; IBorch et all 
I2006t IFaber et alj 120071 : iBrown et all 12007^ Further- 
more, mergers are considered to be the major drivers 
of the growt h of super-massive black holes in the centers 
of ga laxies (jKauffmann & Hachnclt 2000; C roton et al.l 
2006) as they enable large amounts of gas to be funneled 
to, and accrete onto the central BH. 

In the context of the general picture of galaxy evo- 
lution, the following scenario for the evolution of radio 
AGN seems very plausible. High radio power AGN are 
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FlG. 12. — Left panel: The fraction of red galaxies (PI > 0.15; equivalent to U — B > 0.71; left panel), compared to all galaxies, as a 
function of cosmic age for different stellar mass ranges (indicated in the panel). The right panel shows the fraction of radio AGN in the 
red, parent galaxy sample (as defined in the text) as a function of cosmic age for the various stellar mass bins. The fractional change with 
cosmic time has been characterized with a linear slope (a), which is indicated for each stellar mass bin. The radio AGN fraction can be 
related to the average time a galaxy spends as a radio AGN (r ra dio; see text for details), shown on the right-hand y-axis in the right panel. 
Note the different y scales in the left and right panels. Note also that there is no correspondence between the increase of the red galaxy 
fraction, and the decrease of the radio AGN fraction with cosmic time. 



preferably found in galaxies that have not yet reached 
the highest stellar masses and show a strong decline with 
time of their number density. On the other hand, the 
hosts of low power radio AGN are generally old mas- 
sive 'red-and-dead' galaxies and have an almost con- 
stant space density. This implies that the radio trig- 
gering mechanism in a galaxy is a strong function of the 
host galaxy's properties and it is likely that it is linked 
to different stages during the formation of massive old 
galaxies. In this scenario the occurrence of the two radio 
modes would be naturally linked to the overall galaxy 
formation process and their impact and timescales could 
be observationally constrained. 

In particular, the HER AGN phase seems to be closely 
tied to major merger remnants (as many hosts are 
found to be disturbed, have blue colors and show signs 
of on-going star formation; e.g. iHeckman et alj 119861 : 
iBaum et al.l Il992f) . as well as significant black hole 
growth (supported by high - ~ Eddington - accre- 
tion rates observed; see Sec. [5]). These major merg- 
ers are critical for the formation of massive ellipticals. 
On the other hand, the LERAGN phase is mainly ob- 
served in the m ost massive galaxies, where a mecha- 
nism is required (|Croton et al.ll2006t iBower et al.l 12006b 
which prevents intra-cluster gas from cooling onto the 
massive host galaxy. The properties of LERAGN such 
as lo wer amounts of gas and dust compared to HER- 
AGN tide Koff et al |[l99l iLeon et alj[200ll iMiiller et all 
120041; [Evans et alj l2006t lEmonts et all l2008f). a mod- 
est BH gro wth (lEvans et alJl2006t iHopkins et al.ll2006aL 
ICaol 120071 : lHardcastle et al.1 l2007t ). and fully built up 



stellar and BH mass es already at z ~ 1 (Sec. [5l 
iKauffmann et al1l2008l ). support this interpretation. 

In the context of the generally adopted red galaxy 
formation scenario, our results thus suggest that high- 
excitation radio AGN (i.e. high radio-power AGN; 
Li.4GHz > 10 25 W Hz -1 ) occur in a stage of galaxy evo- 
lution presumably after a gas-rich merger which induces 
a phase of significant black hole growth. This active ra- 
dio phase may last 10 6 — 10 7 yr (e.g. Alexander & Leahy 
1987, Shabala et al. 2008). Within the post-merging pro- 
cess the accumulated cold gas will be consumed by the 
host galaxy via star- formation and/or black-hole accre- 
tion. Low-excitation radio AGN (i.e. low radio power 
AGN; Li. 4G hz < 10 25 W Hz -1 ) then represent a later 
stage of the galaxy formation process, when the amount 
of cold gas in the host has already been reduced. In 
this phase the galaxy's BH is accreting quiescently in a 
radiatively inefficient way, i.e. below a certain critical ac- 
cretion rate in Eddington units, a nd most of the released 
energy is in kinetic form (see e.g. lHardcastle et al.l[2007t 
I Merloni fc Heinzl [2008) . This radio phase lasts a few 
times 10 7 yr to a few tim es 10 8 yr fe.g. lMcNamara et al.l 
120051 iNulsen et al.l[2005h . Thus, our results imply that 
high radio power AGN appear in an early stage, while 
low radio power occur in a late stage in the evolution of 
massive galaxies. 

The evolution of the integrated luminosity density, 
shown in Fig. 01 can be well explained in our scenario. 
This diagram can be interpreted as the evolution of 
^i.4GHz for two radio luminous, but fundamentally dif- 
ferent, types of galaxies where the HERAGN and LER- 
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Fig. 13. — Evolution of the comoving 20 cm integrated lumi- 
nosity density (f2i.4GHz) for all radio AGN (filled red curve) and 
star forming galaxies (vertically blue hatched curve) since z = 2.5. 
The evolution of Hi ,4GHz f° r the total AGN population is a super- 
position of the evolution of low ( VLA-COSMOS; ora nge hatched 
curve) and high (model taken from Willott ct al. 2001; red hatched 
region) radio power AGN. The light-gray shaded area of the plot 
shows the redshift range where the evolution of the VLA-COSMOS 
AGN and star forming galaxies has been extrapolated [z > 1.3; see 
text for details). 

AGN reflect modes of vigorous and modest BH growth, 
respectively. 

The luminosity density of the HERAGN equivalent to 
a mode of vigorous BH growth evolves dramatically since 
z = 1.3, consistent with e.g. the evolution of optically and 
X-ray selected quasars (see Willott et al. 2001; see also 
Boyle et al. 2000 and Ueda et al. 2003). The state of mod- 
est black hole growth (in our case the LERAGN) evolves 
significantly slower, and it dominates the total comoving 
AGN radio power output at low redshifts (z < 0.7). At 
these redshifts the most massive, 'red-and-dead', often 
central cluster galaxi es associated with radi o sources, are 
already in place (e.g. iBlakeslee et al.l [20031) thus provid- 
ing an optimal environment for a self-regulating process 
of modest BH growth (at sub-Eddington accretion rates) 
and low radio power AGN activity. 

In summary, the evolution of a radio selected AGN pop- 
ulation seems to be consistent with the possibly bimodal 
evolution of black hole accretion modes (modest vs. vig- 
orous ) that occur during different stages in the process of 
the formation of (massive) galaxies. 

7.2. The co-evolution of star formation and radio AGN 
activity in radio luminous galaxies 

Black hole growth, especially merger driven 
gro wth, is expected to be l i nked to starburst activ- 
ity (ISanders fc Mirabell 119961: IBoyle k Terlevicbl 119981: 
Franceschini et alj Il999t ISandersI l2003t ICroton et all 
20061) . Thus, in Fig. fTS] we compare the evolution of 
the integrated comoving 1.4 GHz luminosity density 
for the VLA-COSMOS star forming galaxies (see S08b 
for details) and the total radio AGN population, likely 



reflecting two separate modes of black hole accretion, 
as argued above. For the latter, f2i.4GHz has been 
obtained by co-adding the contributions of the low- 
power (VLA-COSM OS) and the high radio-power AGN 
(jWillott et al1[2T)0l1 : see Fig. 0). 

As the peak of the cos mic star format ion and AGN 
activity occurs at z > 1 ([Hopkins! l2004h . we extrapo- 
late Oi.4ghz to z = 2.5 as fo l lows. For the high-power 
AGN we use the I Willott et all (|2001l ) model which is well 
constrained b y 3CRR, 6CE, and 7CRS data at z < 3 
(see Fig. 1 in lWillott et"afll2001l) . As can be seen from 
Fig. El the full (weak plus powerful radio AGN) Willott 
et al. model describes well also the evolution of weak ra- 
dio AGN, at least at z < 1.3. Here we assume that this is 
the case also for z > 1.3. Thus, in order to be consistent 
with this model, in which the low-power AGN popula- 
tion is assumed to stop evolving beyond a certain red- 
shift, we assume that the comoving volume densities of 
VLA-COSMOS AGN remains constant for z>1.3. Note 
however that this may not necessarily be realistic if the 
comoving volume density of low- luminosity AGN exhibits 
a turn-over, but at lower redshift compare d to powerful 
radio sources (e.g. Waddin gton et al.1l2001h . 

To extrapolate ^i.4GHz for the star forming galaxies we 
follow Hopkins (2004) who showed, based on a compila- 
tion of numerous studies based on independent data sets, 
that ceasing the evolution beyond z = 2 (given the evolv- 
ing radio luminosity function for star forming galaxies) 
reproduces the observed cosmic star formation history 
well (see e.g. Fig. 1 in Hopkins 2004). Thus, our extrap- 
olation of r2i.4GHz beyond z — 1.3 is fairly robust as it 
is based on independent results drawn from various data 
sets that constrain these cosmic times well. 

The integrated comoving 20 cm luminosity density for 
radio luminous star forming and AGN galaxies, shown 
in Fig. [131 appears to evolve coevally at high redshifts 
(0.7 < z < 2.5), where powerful AGN significantly con- 
tribute to the r2i.4GHz- Both seem to flatten in the 
redshift range of about 1.5 — 2. Although not conclu- 
sive, this is suggestive of a link between the process 
of star formation and radio AGN activity (see al so e.g. 
IBoyle fc Terlevichl Il998t iFranceschini et alj Il999i . This 
can be understood in the scenario of galaxy formation 
where gas-rich mergers of spiral galaxies govern the for- 
mation of gas-poor elliptical galaxies. As the galaxies 
merge, the bulk of their cold gas, that was originally dis- 
tributed throughout the disk of the merging constituents, 
is funneled into the inner few kpc of the merging system 
and fuels both st ar formation a nd AGN activity (pos- 
sibly sequentially: ISandersI 12003). Thus, co-evolution of 
star forming and powerful radio AGN activity is predom- 
inantly expected in systems with significant amounts of 
cold gas, consistent with the properties found in high- 
excitation radio AGN. 

On the other hand, the evolution of the integrated co- 
moving 20 cm luminosity density for star forming and ra- 
dio AGN galaxies decouples below z ~ 0.7 (see Fig. [T3| 
where low radio power AGN with modest BH growth 
dominate the total AGN radio-power output in the uni- 
verse (see Fig. [4]). In the context of our high-to-low- 
radio-power transition scenario proposed in the previous 
section, this AGN activity is triggered in a later evo- 
lutionary stage of a massive galaxy (corresponding to 
LERAGN), where a large fraction of the gas reservoir 
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has already been consumed. Thus, a close link between 
star formation and AGN activity is not expected in this 
regime. Such an expectation is also supported by our 
results on the quenching of star formation in the context 
of the build-up of the red galaxy sequence, derived in 
Sec. 16.31 We have shown that processes causing the star 
formation quenching are not likely to be also governing 
the triggering of low radio power AGN phases. If they 
arc responsible for the decline in the star formation rate 
density since z ~ 1 (at least at the high stellar mass 
end; see e.g. S08b), then no connection between f2i.4GHz 
for low-power radio AGN and star forming galaxies is 
expected. 

We would like to stress that the radio energy injection 
into the surrounding medium of such quiescent modes of 
black hole accretion ha s been invoked in numerous cos- 
mological models (e.g. Croton et al. 2006; B ower et al.l 
I2006t ISiiacki fc Springe]||2006HSiiacki et alj|2007l) as the 
main heating mechanism responsible for preventing fur- 
ther star formation in the given host galaxy (assuming 
that every massive red galaxy undergoes a phase of ra- 
dio AGN activity). This process was the key ingredient 
that allowed a good reproduction of e.g. the high mass 
end of the galaxy stellar mass function in the models. 
To shed light on the plausibility of such a process we 
explore in the next section the injection energies of our 
VLA-COSMOS AGN, and compare them to the recipes 
used in cosmological simulations. 

7.3. The evolution of radio AGN feedback in massive 

galaxies 

In this section we derive the mechanical energy in- 
jected into the surrounding medium by low-power (VLA- 
COSMOS) radio AGN, which are possible candidates 
for preventing star formation in mas sive galaxies as pro- 
posed in recent cosmological models (ICroton et alj|2006t 



Bower et al.l 120061: ISiiacki fc Spri ngcl 2006: Siiack i et al 
20071V and defined as 'radio mode' (see lCroton et al.l2006l 



for details). 

Monochromatic radio luminosity is not a direct indi- 
cator of the mechanical energy output of a radio galaxy 
into its surrounding medium. The latter can be esti- 
mated making use of the observed interplay in galaxy 
clusters between radio galaxies and the hot X-ray emit- 
ting intra-cluster gas, where radio jets induce cavities 
in the hot gas and inflate buoyantly rising bubbles (see 
e-g. iBirzan et alj|2004 I Allen et all 120061: iBirzan et all 
2008). As the total energy within a bubble is the sum of 
its internal energy and the PdV work done by its infla- 
tion, the mechanical radio luminosity is given by the ratio 
of this total energy and the age of the bubble (see e.g. 
iBzrzan et alj|2004 1 Allen et alj|2006t iBzrzan et alll2008f). 
Using mechanical luminosities derived by iBirzan et al.l 
(2004) for a sample of 18 radio galaxies in clusters or 
groups, Best et al. (2006) have found that the relation 
of mechanical and monochromatic radio power is best 
described by a power-law, however with a large scatter 
(see Fig. 1 and eq. 2 in Best et al. 2004). Based on this 
relation, they have computed the total comoving volume 
averaged mechanical heating rate, ^Lmoch, in the local 
universe provided by radio luminous AGN, which was 
found to be a factor of 10-20 lower than the local den- 
sity of radio mode heating predicted in the Croton et al. 



(2006) semi-analytic model. 

Here we follow the prescription of lBest et al.l (|2006f ) to 
estimate the evolution of OLmechj using the most recent 
calibration for the c onversion of mono chromatic radio to 
mechanical power (iBzrzan et al.l f2008h . First, at each 
redshift we convert the radio AGN monochromatic lu- 
minosity function (given by our best fit PLE and PDE 
models; see Sec. 14.11) into a mechanical luminosity den- 
sity function using eq. 16, given bv lBzrzan etall ((2008). 
Then we integrate it above a minimum mechanical lu- 
minosity equivalent to Li. 4 gh z = 10 21 W Hz -1 . The 
resulting evolution of f^Lmcch is shown in Fig. 1 141 where 
we compare it to the density of the radio mode heat- 
ing given in the Croton et al. (2006) cosmological model. 
The latter has been obtained from Fig. 3 in Croton et al., 
where we converted their comoving black hole accretion 
rate density, meH (shown on the right-hand side of the 
y axis), into the volume- averaged heating rate utilizing 
fiLmcch = ^™bh c 2 ; where r\ is the standard efficiency of 
the conversion of mass into ener gy (77 = 0.1), and c i s the 
speed of light (see also eq. 11 of lCroton et al.ll2~006f) . 
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Fig. 14. — The cosmic evolution of the volume averaged me- 
chanical heating rate flLmech for low-power VLA-COSMOS AGN, 
which arc likely candidates for the 'radio mode' heating invoked 
in cosmological models (filled curve). The volume averaged ac- 
cretion rate is shown on the right-hand side y-axis (see text for 
details). The uncertainties in the correlation between Li.4qhz and 
Lmcch (given inl Birzan et al. | 206g) are illustrated by the light-gray 
shaded area. Also shown is the evolution predicted in the Croton 
et al. (2006) semi-analytic model (thick line), and its possible value 
lowered by a factor of 2-3 (indicated by arrows) due to a system- 
atic over-estimation of the heating rate in the model (see text for 
details). 

It is important to note that both the model- and 
observation-based heating rates are subject to large un- 
certainties. The model cooling rates are likely to be over- 
estimated by a factor of 2-3 (indicated in Fig. [T4")) due to 
the nature of the underlying cooling flow model (Best et 
al. 2006). On the other hand, the observationally derived 
heating rates have uncertainties of the order of ~ 0.8 dex 
due to the scatter in the conversion of monochromatic to 
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mechanical radio luminosities. Furthermore, they may 
be strongly under-estimated (up to a factor of 6) as the 
work done during the inflation of a bubble may be signif- 
icantly higher than the reversible pdV work usually as- 
sumed in the computation of mechanical luminosities (see 
Binney et al. 2007 for details). In addition, the minimum 
mechanical power equivalent to Li.4ghz = 10 21 W Hz -1 
has been chosen rather arbitrarily. For example, a lower 
threshold would yield a higher ^Lmcch- Nonetheless, the 
qualitative agreement between the model and observa- 
tions seen in Fig. [14] is encouraging for the idea that ra- 
dio luminous AGN play an important role in the process 
of galaxy evolution. 

Based on two samples of local massive red galaxies 
with avail a ble o ptical, radio and X-ray observations, 
iBest et all (|2006t ) have shown that radio source activity 
may indeed control the growth rate of the host galaxy, 
as has been proposed in cosmological models. Best et al. 
have found that, across a wide range of stellar masses, 
the radio source heating rate agrees remarkably well with 
the cooling rate of the host galaxy's surrounding hot gas 
halo. Averaged over time, the recurring radio source ac- 
tivity can then balance the radiative losses from its sur- 
rounding hot gas halo and thereby suppress the cooling 
of the gas, and in this way control the growth rate of the 
galaxy (see e.g. lCroton et al . 2006). Here we have shown 
that various properties of z ~ 1 (low power) radio AGN 
are comparable to the properties of local radio AGN. 
Thus it is very likely that the same radio heating - hot 
gas cooling interplay, observed locally, takes place also in 
z ~ 1 radio AGN. This then strongly argues in favor of 
the plausibility that (low-power) radio AGN activity (i.e. 
'radio mode') is a key factor that prevents further mass 
growth in massive 'red-and-dead' galaxies across cosmic 
times, at least since ~ 4.5 Gyr after the Big Bang. 

8. SUMMARY AND CONCLUSIONS 

Using the largest sample of low luminosity radio AGN 
at intermediate redshifts available to-date drawn from 
the VLA-COSMOS survey, we explored the evolution 
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and composition of the radio AGN population out to a 
rcdshift of z~l. In particular, we derived the radio AGN 
luminosity function in several rcdshift bins. Compared 
to powerful (Li. 4GH z > 5 x 10 25 W Hz" 1 ) radio AGN 
the low radio power (Li. 4GH z < 5 x 10 25 W Hz" 1 ) AGN 
probed here show only a weak decline in their luminosity 
density since z = 1.3. 

An analysis of the host properties of weak and powerful 
radio AGN shows that they form two relatively distinct 
classes, with powerful radio AGN being related to modes 
of vigorous BH growth whereas weak radio AGN are re- 
lated to only modest BH growth. This is also reflected 
in the types of host galaxies with weak radio AGN be- 
ing preferentially found in the most massive and evolved 
galaxies (even out to z^l). 

Joint analysis of the evolution of the 20 cm luminosity 
density of the AGN and the star formation rate density 
suggests that the powerful radio sources are preferably 
triggered in major mergers following shortly after the 
phase of massive star formation at high redshifts. This 
would also naturally explain the strong decline of these 
sources below a redshift of z=l. As the volume density 
of weak radio AGN stays fairly constant, they can signifi- 
cantly contribute to the heating of the their surrounding 
medium and thus inhibit gas accretion onto their host 
galaxies as recently suggested for the radio mode in cos- 
mological models. 
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